A 3D topological insulator quantum dot 
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3D topological insulators (TIs) are materials with topologically protected interface/surface states 
of massless Weyl fermions. This protection is manifested by the suppression of the backscattering 
caused by nonmagnetic impurities and edges on the surfaces^^ and by the robustness against any 
type of surface modifications. Consequently, 3D TI nanostructures are very interesting because of 
their large surface-to-volume ratio. The topological surface states in nanostructures exhibit a phase 
coherence length of several hundred nanometer d 11 ! 12 ! Experiments on both the physical and chemical 
synthesis of TI nanostructures were performed recently to understand their transport properties at 
the nanoscale^^. Here we show the model of a quantum dot (QD) consisting of a spherical core- 
bulk heterostructure made of 3D TI materials, such as PbTe/Pbo.3iSno.69Te, with bound massless 
and helical Weyl states existing at the interface and being confined in all three dimensions. The 
number of bound states can be controlled by tuning the size of the QD and the magnitude of the 
core and bulk energy gaps, which determine the confining potential. We demonstrate that such 
bound Weyl states can be realized for QD sizes of few nanometers. We identify the spin locking 
and the Kramers pairs, both hallmarks of 3D TIs. In contrast to topologically trivial semiconductor 
QDs, the confined massless Weyl states in 3D TI QDs are localized at the interface of the QD 
and ex hibit a mirror symmetry in the energy spectrum. Because of the possibility to generate spin 
current^™ the 3D TI QD opens up avenues to develop nanostructures for spintronics devices. 
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3D topological insulators (TIs) are narrow-bandgap 
materials with topologically protected gapless (metallic) 
surface / inte rfac e states that appear within the bulk in- 
sulating gap^l These states are characterized by the 
linear excitation energy of massless Weyl fermions. In 
such materials, the spins of the Kramers pairs are locked 
at a right angle to their momenta on the Fermi surface 
due to spin-orbit coupling. The surface states are pro- 
tected by time reversal symmetry, making them robust 
against perturbation and exhibiting the characteristic 
suppression of backscattering from edges and nonmag- 
netic impuritieiP^. Unlike in normal metals with high 
disorder, where scattering is so strong that electrons get 
localized due to Anderson localization, on a topologically 
protected surface Anderson localization and backscatter- 
ing are suppressed, resulting in ballistic and dissipation- 
less charge transport. Such states are of great importance 
in low-power opto-spintronics that use spin currents un- 
der efficient electrical control 5 -^. Decoherence can be cir- 
cumvented by highly polarized spin states with helical 
spin texture^QSl The presence of the strong spin-orbit 
coupling on the surface of 3D TIs makes these materials 
excellent spin current generators. 

The topological surface states in nanostructures 
have a p hase coherence length of several hundred 
nanometer d 11 ! 12 !. Studies have been done to explore 
them as poss ible ca ndidates for spin current generators 
in spintronic a 6 l 13 l 14 [ Suppression of backscattering be- 
tween states of opposite momentum and opposite spin 
was clearly observed in an experiment on a random al- 
loy Bii-^Sbj; with strong atomic scale disorders. This 
observation further confirms the potential use of 3D TIs 
for coherent spin transport phenomena. 

The electronic properties of the surface states of 3D 



TIs can be analyzed using Angle Resolved Photoemis- 
sion Spectroscopy (ARPES)^£"H] anc i scanning tunnel- 
ing spectroscop y^ 5 * 1 ^ Owing to their technological ap- 
plications as well as fundamental concepts in solid state 
physics, these surface states, despite co ntribu tions com- 
ing from the conducting bulk band stateiP^R are probed 
and studied in the spin transport regime 22 . The sur- 
face states are probed in an experiment where the size 
of the Fermi surface measured by Shubnikov-de Haas os- 
cillations are matched to the carrier density obtained in 
low- field Hall measurement!^. In such experiments, ap- 
plication of a moderate to strong magnetic field reduces 
the conducting bulk bands' contribution to the conduc- 
tivity. These attempts have been focused to get pure spin 
currents that have origin on the metallic surface states. 
To measure spin currents, TIs can be driven out of equi- 
librium by illuminating them with a circularly polarized 
light. The direction of the resulting photocurrents car- 
ried by helical Weyl fermions can be altered by changing 
the helicity of light from left circularly polarized to right 
circularly polarized light and vice versa^^ES 

An important feature of topological surfaces states is 
that they are robust against any type of surface modifica- 
tions, in contrast to topologically trivial semiconductors 
where dangling bonds and reconstructions are common 
on the surfaces and modify largely their surface electronic 
properties. In 3D TI nanostructures the special proper- 
ties of topologically protected surface states of TIs are 
very interesting because the large surface-to- volume ratio 
amplifies any physical effect originating from the surface. 
In such structures, the contribution to the conductivity 
from surface carriers is much greater than that from bulk 
crystals. In addition, the chemical potential can be elec- 
trically tuned using a gate voltage. For example, in an 
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interference typ e exp eriment such as Aharonov-Bohm 
(AB) oscillation d 25 * 26 !, the coherent propagation of the 
Weyl electrons around the perimeter of a nanoribbon pro- 
vides excellent evidence of the topological nature of the 
surface states in TI nanostructures^. Experiments on 
both the physical and chemical synthesis of TI nanostruc- 
tures have been done recently to understand their trans- 
port properties at the nanoscale^HSD Recently, in a TI 
QD with tunable barriers based on ultrathin Bi2Se3 films. 
Coulomb blockade with >5 meV charging energy was ob- 
served and evidence of tunneling into excited states was 
providecP4 

So far, a theoretical study of electronic properties of 
2D helical states occurring at the nanoscale of 3D TIs, 
such as in QDs, is still lacking. In this article, we present 
for the first time the study of bound Weyl states that 
are confined at the interface of a spherical core-bulk 
heterostructure QD made of 3D TI materials such as 
Pbi-xSn^Te. We show that at the interface massless 
Weyl fermions are confined in all three dimensions. The 
directions of spin and momentum are tangent to the sur- 
face of the QD, as shown in Fig. [3] Remarkably, their 
inherent spin-momentum locking property exists even in 
the three-dimensional confinement of a QD. Since the 
massless Weyl fermions have a linear dispersion, there is 
a mirror symmetry in the energy spectrum between pos- 
itive and negative energy states, in contrast to topologi- 
cally trivial semiconductors where holes in p-type valence 
bands are typically heavier than electrons in s-type con- 
duction bands, resulting in completely different energy 
spectra. We demonstrate that this symmetry in energy 
spectrum is preserved for the QD spectrum. 

In Fig. [T| we show the model of our spherically sym- 
metric 3D TI QD of a core-bulk structure with a single 
interface at radius r = r$ [see Fig. [l]a]. This core-bulk 
structure consists, for example, of an inner core of PbTe 
and an outer bulk of Pbo.31Sno.69Te with bandgaps of 
0.187 and -0.187 eV, respectively, or vice versa, so that 
Weyl fermions are generated at the interface^. Here we 
used the bandgap formula provided in Ref. l3"2Tbr de- 
termining x. Note that the band crossing happens in 
Pbi^Sn^Te at x = 0.35 at 4 K. The Weyl fermions 
are subjected to the spherically symmetric potential A(r) 
(Fig. [l]b), which breaks the translational invariance in 
radial direction. 

To understand the energy spectrum and helical prop- 
erties of massless Weyl fermions at the interface of a QD , 
we start with a Dirac Hamiltonian within the envelope 
function approximation. Neglecting the far band terms, 
we have 



H 



A(r) 
cr.p 



cr.p 

-A(r) 



(1) 



where cr are the Pauli matrices, p is the momentum oper- 
ator with Fermi velocities t^and v\i in angular and radial 
direction, respectively, A(r) = e g (r) /2 is the gap en- 
ergy parameter. It is important to note that the spheri- 





Figure 1: A heterostructure spherical core-bulk 3D TI QD 
with a single interface, a. The 2D interface of the 3D TIs 
can be curved to make a spherical core-bulk QD. Inner and 
outer materials are chosen in an appropriate ratio so that 
Weyl fermions are generated at the interface. The arrows 
at the corners are drawn to show the infinite size of the 
host. The core and bulk host can be chosen as PbTe and 
Pbo.31Sno.69Te or vice versa, b. The potential A(r) with 
A (r — ro) = —A (ro — r) breaks the translational invariance 
of the crystal in radial direction. Weyl fermions possess zero 
mass; therefore, 3D TI QDs exhibit a mirror symmetry be- 
tween positive and negative energy states. The energy of the 
bound interface states depends on the size of the QD and the 
strength of the potential. As an example, two bound states 
at the interface are shown with energies +e and — e (short 
dashed lines) for a QD of size ro =2 nm. 



cally symmetric A (r) breaks the translational invariance 
of the crystal in radial direction and has the symmetry 
A (r — ro) = —A (ro — r), where ro is the radius of the 
QD. The eigenfunctions of H are four-component spinors 

X- 



where X- an d x+ are two-component spinors 



corresponding to the L_ and i + band, respectively, such 
as in Pbi-^Sn^Te. There are two important points: first, 
we assume our QD has spherical symmetry. This means 
the angular parts are separated from the radial part of 
the Dirac Hamiltonian given by equation 0. Second, in 
place of a constant mass term in the Dirac Hamiltonian, 
we have a spherically symmetric potential term A (r) . 
Since this system has spherical symmetry, we follow the 
derivation of the solution for the central-force problem 
of a hydrogen atom in relativistic quantum mechanics, 
as treated explicitly in Ref. 33. After eliminating the 
angular parts, the radial part of the Dirac Hamiltonian 
given in equation ([I]) takes the form 



H = 



A(r) 



-A(r) 



(2) 



where vu = 2.24 x 10 5 m/s for Pbi^^Sn^Te and k — 
± (j + J;) is a nonzero integer which can be positive or 
negative, j being the total angular momentum quan- 
tum number. For given k, it is known from relativis- 
tic quantum mechanics that the angular momenta l_ 
and Z+ for x- an d X+ are determined by the relations 
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-« = 3 U + (l- + 1) + 1/4 and K = j (j + 1) - 

(Z+ + 1) + 1/4, respectively. Now we solve the eigen- 
value equation by squaring the Hamiltonian given by 
equation pj, resulting in the new eigenvalue equation 
H 2 ip — e 2 ^, which consists of the two coupled second- 
order differential equations 



dr 2 



(3) 

where /± = r X ±, P = l/v\\h and A = /V( A o - e 2 )- It 
is to be noted that the parameter A behaves like a wave 
vector k. whose allowed quantized values determine the 
particle's energy levels. The set of equations (|3| are cou- 
pled inhomogeneous differential equations which can be 
solved by the standard Green's function technique. In a 
flat geometry such as the thin layer of a 3D TI, the lin- 
ear Weyl spectrum k± = is approximately equal to the 
soliton spectrum in the ID Peierl's insulator and A (z) 
can be chosen to be A(z) — A (oo) tanh (z/l)™The cor- 
responding interface states are localized along the z-axis. 
Consequently, we adopt a similar potential along the ra- 
dial direction of the form A (r ) = A sgn (r — r^j , 

where A D is the strength of the potential and r is 
the radius of the spherically symmetric interface of the 
QD. Hence, the source term in equation ^ is simply 

F(r ) = 2A j3f±(r )r 2 5 (r — r D ^j . The set of equations 

^ can be solved by using the corresponding differential 
equation for the Green's function, i.e. 



'/ (V-;/ 1 

dr \ dr 



V 



k{k± 1)) 



G T =6(r- r 



The solutions regular at r = with outgoing wave be- 
havior at r — > oo are the product of spherical modified 
Bessel functions of the order k for G_ and of the order 

k — 1 for G+, i.e. G_ (r, r , a) = CJI K ( Ar< ) K K ( Ar> ) , 
G+ (r, r , a) = C + l n -i (Ar<)/C re _i (Ar>), where r< 

(r>) is the smaller (larger) of r and r . The func- 
tions I{Xr) and K, (Ar) are, respectively, the first and 
the second kind of modified spherical Bessel functions, 
and Czf are the normalization constants. These con- 
stants are determined by the discontinuity in slope im- 
plied by the delta function in equation Q. Integration 
is performed at the interface of the QD along the ra- 

r +n 

1, where r\ is an infinites- 



dial direction: 



dr 



v, 



r. Con- 



imal quantity with T] > 0. For r = 
r < = r and for r = r — r/, r> = r . 
sequently, the normalization constants are determined to 
be: C_ = l/Xr'oW K and G+ = l/Ar^W«_i, where = 
T« (Ar') < (Ar) (Ar)/C K (Ar')] and W K -i = 

l K -t (Xr) <_! (Ar) (Ar)/C K _! (Xr 



the Wronskians of I(Ar) and IC(Xr), respectively, for 
k and k—1 order, and X (Ar) and /C (Ar) are deriva- 
tives of the Bessel functions. The Wronskian of two 
linearly independent functions is proportional to 1/r 2 
for Sturm-Liouville type equations such as equation Q 
(see Method A). The solutions of the inhomogeneous 

equations (3) are = J Gzp (r, r , Xj F(r )dr — 

2A„/3 J G T 7r, r , A^j /±r 5 (r — r ^j dr . Using the so- 
lutions for the Green's functions and integrating we ob- 
tain 



/_ = 2A /3/ + I K (Ar<) K H (Ar>) / XW R 



(5) 



/+ = 2A /3/_I K _! (Ar<)/C K _! (Ar>) /AW K _! (6) 

where r < (r>) is now the smaller (larger) of r and r . A 
transcendental equation is obtained by solving equations 
^ and (|6| and evaluating at r = ro, 



[zZ K (z) (z)] [zZ K _! (z) JC K -! (z)] = l/4A 2 /3 2 r 2 (7) 

where z = Ar . Each term in the square bracket on the 
left hand side of equation ^ is a monotonically decreas- 
ing function of z (for z > 0), with maximum value of 
1/(2k + 1) for K th order term and 1 / (2k — 1) for (k — l) th 
order term occurring at z = (see Method B). There- 
fore, their product has a maximum value of 1/(4k 2 — 1) 
at z = and is equal to l/4A 2 /3 2 r 2 . For z — > oo, 
[zI K (z) fC K (z]j — > l/2z. In Fig. [2] we show the plot 
of equation (l7| where the function F (z) is defined as 
F(z) = [zl K (z)lC K (z)] [zZ K _! (z)/C K _! (z)]. Since F (z) 
is a monotonically decreasing function, for each k. there is 
at most a single solution given by the intersection of F (z) 
with the constant l/4A 2 /3 2 r 2 (dashed line and solid line 
in Fig. |2j. The critical limit for having a single solution 
is determined by the intersection at the maximum value 
of F (z), which occurs at z = 0. This means that there 
exists a single solution of equation |7]) for each k as long 
as the condition l/4A 2 /3 2 r 2 < 1/(4k 2 — 1) is satisfied. 
Fig. [2] shows the plot of the first three different values of 
k, k = 1 (red), 2 (blue) and 3 (pink), each a monotoni- 
cally decreasing line (solid line) cut by a horizontal line 
(dashed line) at most one time. Since A = ftyj (A 2 , — e 2 ), 
each single solution gives rise to two bound states with 
same magnitude but opposite sign of energy. Hence, if 
the Fermi energy is set at the origin where the bands 
cross, there is an electron-hole mirror symmetry in the 
energy spectrum. Indeed, this makes sense since Weyl 
fcrmions are massless at zero band gap with the linear 
dispersion relation. This result is in stark contrast to 
topologically trivial semiconductor QDs where holes are 
typically heavier than electrons and therefore break the 
electron-hole mirror symmetry in the energy spectrum. 

As the size of the QD grows, it is filled with more and 
more bound states (see Fig. [2]) where for smaller value 
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Figure 2: Plot of equation Q. The function F (z) is defined as 
F (z) — [zX K (z) K.k{z)][zZ k -i (z) K. K -i (z)]. Intersection of 
horizontal (black) dashed line and monotonically decreasing 
solid line at z = gives the minimum threshold of the size of 
a QD to have two bound states, one positive and one negative 
energy state, for a given confining potential. For the larger 
size of the QD (smaller value of F (z)), multiple bound states 
are available, corresponding to multiple intersection points. 
As the size of the QD increases, the available positions of these 
bound states lie deep inside the confining potential well. The 
intersection points A, B and C are example points where we 
evaluate the wavefunctions. The energy of the bound states 
are determined by the relation z = \r . 
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Figure 3: This schematic shows the partial momenta p_ and 
p+ of the X- an d X+ components on the surface of the QD 
which give rise to the angular momenta 1_ and 1+. In this 
visualization it is obvious that the spin s is locked perpen- 
dicular to the momenta p_ and p+ and perpendicular to the 
interface, which is a hallmark of 3D TIs. 
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Table I: X- an d X+ components 



of F (z). a horizontal dashed line makes multiple cuts at 
different values of the energy (i.e. z) for different n. It is 
important to note that for negative k, the solutions di- 
verge at the origin and are therefore physically not valid. 
This result has profound implications because the sign 
of k determines whether j is parallel or antiparallel to 
the spin s (see Ref. I3"3")l . Since k is only allowed to be 
positive, only one spin orientation with respect to j is 
permitted. This corresponds to the spin locking effect, 
which is a hallmark of 3D TIs (see Fig. |3| . The condition 
1/4Aq/3 2 t-q = 1/(4k 2 — 1) determines the lower limit of 
the size of the QD to hold two bound interface states, 
a positive and a negative energy state, for a given value 
of the confining potential strength. The critical QD size 
depends on the Fermi velocities and band gaps of the 3D 
TI materials. In Pbi-^Sn^Te, the strength A G of the 
confining potential can be taken as half of the band gap, 
by definition, of PbTe, which is 0.0935 eV. Taking the 
Fermi velocity on the surface to be 2.24 x 10 5 m/sJ231 re- 
sults in a critical QD size of ro = 1.4 nm for k = 1 at 
z = 0. Similarly for k — 2 at z = 0, the critical QD size 
for Pbi_a.Sn.rTe is ro = 3 nm. The energy of the bound 
states are determined from z = Ar D , which gives a very 
shallow energy level of e = ±A Q for z = 0. 

For a given value of k, quantum numbers character- 



izing the wavefunctions X- an d X+ can be determined. 
For k = 1,2,3 and 4, the possible combination of 
the quantum numbers are shown in Table [I] for both 
spinors x- and \+- Note that only the positive values 
of k are shown since a negative k is prohibited. Here 
we observe that the x~ component is characterized by 
the spin being antiparallel to its angular momentum, 
whereas the x+ component is characterized by the spin 
being parallel to its angular momentum. We show 
now how to identify the Kramers pairs. According 
to Kramers theorem, which applies to a time-reversal 
invariant system, a spin 1/2 state is at least twofold 
degenerate on the surface of a 3D TI. Hence, using the 
4-spinor notation [|x— ) j |x+)L spin-orbit kets \mi,m s ), 
and Clebsch-Gordan coefficients, we obtain the following 
examples of Kramers pairs. For k = 1, a 4-spinor state 

0, i)") , |0, |) with m 
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Kramers partner 
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calculated for the intersection point A (ref. Figure 2) 
I 



calculated for the intersection: point B 
(ref. Figure 2) 1 




Figure 4: Spatial dependence of the two-component spinors 
X-and x+ inside and outside the QD calculated for the in- 
tersection point A shown in Fig. [2] The QD is made of 
PbTe/Pbo.3iSno.69Te core-bulk heterostructure with a core 
of size ro = 2 nm. The vertical dashed line denotes the posi- 
tion of the QD boundary. The solid horizontal lines represent 
the energy eigenvalues e± = ±0.8A o , as marked on the right 
vertical axis. 
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respectively. A state with l z 



= 2 



and 

M) 

has 

two pairs of Kramers partner and a state with = 3 
has 3 pairs Kramers partner. In general, the number of 
Kramers pairs is determined by the spin multiplicity for 
each rrij value. 

In Figs. [4] and [5] we show the spatial wavefunc- 
tions of the \- an d X+ components inside and out- 
side the QD made of the core-bulk heterostructure 
PbTe/Pbo.3iSno.6gTe. Figs. |4]shows the example of the 
intersection point A (see Fig. |2j, and Fig. [5] shows the 
example of the intersection points B and C (see Fig. |2| . 
Since the 4-spinors must be continuous at the boundary, 
also each of the 2-spinor components must be continuous, 
i.e. x 1 - an d X+ must be equal to the spinors and 
X+ Ut , respectively, at the boundary of the QD surface. 
The horizontal solid and short dashed lines in Figs. [4] 
and [5] represent the energy eigenvalues, respectively, at 
the intersection point A, corresponding to a size of the 
QD of 2 nm, and at the intersection point B and C, cor- 
responding to a size of the QD of 3.5 nm. Eigenvalues 
are given by e± — ±0.80A o at point A, e± = ±0.9lA o 
at point B, and e± = ±0.48A o at point C. 

In conclusion, we have shown that Weyl fermions can 
be confined in all three dimensions at the spherically 
shaped interface between two narrow-bandgap semicon- 
ductor alloys, such as the core-bulk heterostructure made 
of PbTe/Pbo.3iSno.6gTe. This configuration provides us 



0.30 
2 0.25 

to 
c 

^0.20 



co 

£ 0.15 

0! 



0.10 



0.05 



0.00 



point C 



1 






1 

A\ 

ii'A 








inside / i , 

/&•• 


i outside 




/ i/ Y 1 

1 M \ 


\ \ 
* \ 
% \ 
\ \ 




•1 '/ c 1 

'/ , / ° 
'/ tl 


\ * >«. 




.'/.'/ 1 


\\» 




* f * f 
/ /// 1 









1.0 



0.5 



< 

*4 

o 



0.0: 



-0.5 



r o 4 
r (nm) 



-1.0 



10 



Figure 5: Spatial dependence of the two-component spinors 
X-and x+ inside and outside the QD calculated for the inter- 
section points B and C shown in Fig. [2] The QD is made of 
PbTe/Pbo.3iSno.69Te core-bulk heterostructure with a core 
of size ro = 3.5 nm. The vertical dashed line denotes the 
position of the QD boundary. The solid horizontal lines rep- 
resent the energy eigenvalues e± = ±0.91A o at point B and 
e± — ±0.48A o at point C, as marked on the right vertical 
axis. 



with the model of a spherical 3D TI QD with tunable size 
ro and potential An, which allows for complete control 
over the number of bound interface states. In contrast 
to topologically trivial semiconductor QDs, the energy 
spectrum of the 3D TI QD exhibits a mirror symme- 
try between positive and negative energy states, which 
is typical for massless Weyl fermions. While in trivial 
semiconductor QDs the states are confined inside the core 
material, in a 3D TI QD the states are confined at the in- 
terface between the core and the bulk materials. We have 
shown that each of these bound states are characterized 
by the intrinsic spin-momentum locking property with 
helical spin states, which they inherit from the 2D inter- 
face structure. Therefore, 3D TI QDs have completely 
new physical properties, which might be very useful for 
the development of future spintronics devices due to their 
inherent topological protection against backscattering. 
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Methods 



B. Maximum of the function F (z) 



A. Calculation of the Wronskian 

The Wronskian of the functions I K (z) and K. K (z) is 
defined aiP^ 

W K [1 K (z) , K K (z)} = 1 K (z) K.' K (z) - T K (z) K K (z) , (8) 

where the prime denotes the derivative of the function. 
For independent solutions, it is to be noted that Wron- 
skian is proportional to 1/p (x) in a Sturm-Liouville type 



dx 



p (x) J 2 + g(x)y = 0. Therefore, Wron- 



equation 

skians in the text are calculated to be 



and 



W K [l K (z),JC K (z)} 



W K . 1 [l K . 1 {z) ) ]C K - 1 {z)} = — 



(9) 



(10) 



The limiting forms of modified Bessel functions are 
given by 



Ik ( z ) — r(«+i) (2) 
K K (z) = ^ (§) 



k asz^O. (11) 



The modified spherical Bessel functions can be written 
in terms of modified Bessel functions as 

Therefore, the function 

F (z) = [zT K (z) K K {z)\ \zX K _ x (z) K K _ X {z)\ (13) 
has the limiting form F [z) = . l_ x as z — > 0. 
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